Half-Integer Filling Factor States in Quantum Dots 



O 
O 



A. Harju,^^0 H. Saarikoski,^' ^ and E. Rasanen^^^ 

'Laboratory of Physics, Helsinki University of Technology, P.O. Box 1100, FI-02015 HUT, Finland 

^Kavli Institute of NanoScience, Delft University of Technology, 2628 CJ Delft, the Netherlands 
^Institut fiir Theoretische Physik, Freie Universitdt Berlin, Arnimallee H, D-14195 Berlin, Germany 

(Dated: February 6, 2008) 

Emergence of half-integer filling factor states, such as v—5/2 and 7/2, is found in quantum dots by 
using numerical many-electron methods. These states have interesting similarities and differences 
with their counterstates found in the two-dimensional electron gas. The i/ = 1/2 states in quantum 
dots are shown to have high overlaps with the composite fermion states. The lower overlap of the 
Pfaffian state indicates that electrons might not be paired in quantum dot geometry. The predicted 
u = 5/2 state has high spin polarization which may have impact on the spin transport through 
quantum dot devices. PACS: 71.10.-w, 73.21.La, 73.43.-f, 85.35. Be 
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The fractional quantum Hall state at filling factor 
1/ = 5/2 in the two-dimensional electron gas (2DEG) 
is widely thought to be a Moore-Read state of p-wave 
paired electrons with unusual topological properties 0. 
The state has recently attracted much theoretical inter- 
est together with other half-integer filling factor states, 
mainly due to proposals for using it in quantum comput- 
ing y . The 2DEG i/ = 5/2 state is expected to consist of 
two Landau levels (LL), where the filled lowest LL (OLL) 
is spin-compensated and the half-filled next LL (ILL) is 
spin polarized state described by a Pfaffian (PF) wave 
function This wave function explicitly describes the 
electron pairing which is analogous to that in the BCS 
theory of superconductivity. Due to the half-filling of 
ILL at v = 5/2, the properties of the i/ = 5/2 state is 
usually analyzed using the v = 1/2 state with an inter- 
action Hamiltonian corrected for the screening effects of 
the electrons in the filled OLL. 

Unlike half- integer filling factor states in the 2DEG, 
the corresponding states in quantum dots (QD's) have 
got almost no attention until now (see Ref. 4 for an excep- 
tion). In this Letter we analyze the electronic structure 
of QD's in magnetic fields and find new states which can 
be regarded as finite size counterparts of the half-integer 
filling factor states in the 2DEG. Since the emergence 
of fractional quantum Hall states is always a result of 
electron-electron interactions, we use numerical methods 
which can take into account the complex correlation ef- 
fects in the system. The results indicate analogies but 
also important differences between the states in the half- 
integer filling factor range in the 2DEG and QD's. The 
mean-field method predicts that due to the presence of 
an external confining potential in QD's, the v = 5/2 state 
appears as a mixture of = 2 and 3. We analyze also 
the internal structure of the many-body state at = 1/2 
in QD's and compare it to the PF wave function and 
the composite fermion (CF) model Q. Unlike the CF 
model, which yields high overlaps with the exact state, 
the PF wave function is found to poorly describe the pure 
I' = 1/2 state and also the excited ILL of the ly — 5/2 
state in QD's. Therefore, we find no evidence of electron 



pairing near half-integer LL fillings in QD's. We finally 
discuss the impact of the i/ = 5/2 state on the ground 
state energetics and on the spin transport through the 
QD and find that signatures of this state can be found, 
e.g., in magnetization measurements or electron tunnel- 
ing experiments. 

We define our QD system by an effective-mass Hamil- 
tonian for N electrons 
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with material parameters for GaAs, i.e, m* = 0.067, 
e/eo = 12.4, and g* — 0.44. Unless stated otherwise, we 
apply a parabolic confining potential, Vc{r) — m*ujQr'^ /2 
with a strength fkjQ = 5.7iV^^/^ meV. The N depen- 
dence in ujQ is introduced to keep the electron density in 
the dot approximately constant. We solve the associated 
Schrodinger equation using the mean-field spin-density 
functional theory (SDFT) and the exact diagonalization 
(ED) method 0. We use the SDFT to calculate the elec- 
tronic structure of large quantum dots which are beyond 
reach of exact many-body method. The internal struc- 
ture of the half-filled LL is then analyzed in detail with 
the ED method which is, however, restricted to fairly low 
electron numbers. 

We analyze the polarization and occupations of the 
LL's in large QD's using the SDFT. We focus here on the 
filling factor regime v > 2. One finds that the filling fac- 
tor for the finite electron systems is position dependent 
7]. Even at integer filling factors (where most physics 
can be understood on a single-particle level), the 2DEG 
and QD's have differences due to the overlapping of LL's 
in QD's in low magnetic fields. In QD's we identify the 
integer v as states, where the increasing magnetic field 
has emptied one of the LL's completely of electrons. For 
example, v — 3 corresponds to case where magnetic field 
has emptied 2LL and the spin-compensated ILL has ap- 
proximately one third of the electrons, the rest being at 
OLL. In the 2DEG, ILL is spin split at = 3 and has 
one spin channel fully occupied. 
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FIG. 1; (color online) (a) Ground-state spin of quantum dots 
with A'' — 12, 30, 48, and 60 electrons, respectively. The 
magnetic fields are scaled linearly to filling factor range 2 < 
1/ < 3. (b) Spin densities of the A'' = 60 quantum dot at v — 
5/2. The spin-polarized second Landau level gives rise to high 
spin polarization in the core region of the dot. (c) Occupation 
of the single-electron (Kohn-Sham) states on different Landau 
levels in the A'^ = 60 quantum dot at2<i^<7/2. 



Figurenja) shows the calculated spin polarization S = 
(jVT - N^)/2 of the ground state of QD's with N = 12, 
30, 48, and 60 electrons, respectively. Our starting point 
is the v = 3 state discussed above. When the magnetic 
field is increased electrons move from ILL to OLL. This 
occurs mainly for spins anti-parallel to the magnetic field. 
We find a peak in S between v = 2 and v = 3 where all 
the spins anti-parallel to the magnetic field have fallen 
to OLL leaving ILL totally spin polarized with angular 
momentum states / = — 1, 0, 1, . . . , will — 2. Since the 
spin polarization of ILL agrees with the expected spin 
polarization of the Moore- Read state in the 2DEG, we 
identify this state as a finite-size counterpart oi v — 5/2. 
In addition, i/avg = 2.54 (see Ref. 0) for this state when 
N = 60. The spin density of v = 5/2 state at iV = 60 
shows that the spin polarization is concentrated at the 
core of the dot [Fig. H^b)]. This result is in contrast 
with the structure of the 2DEG where the single-particle 
states on ILL are uniformly occupied, whereas in QD's 
the occupation has a compact structure due to the exter- 
nal confinement. Thus, the ly — 5/2 state in QD's can be 
interpreted to be composed of two different filling factor 
domains, i.e., v = 2 ai the edge and z/ = 3 at the center. 
We find also an analogous state in lower magnetic fields 
1/ > 3 where the 2LL is spin polarized. We identify this 
state a.s v = 7/2. The Kohn-Sham occupations of LL's 
for 1/ = 2, 5/2, 3, and 7/2, respectively, show the mech- 



anism of LL filling in QD's [Fig. ^c)]. Due to partial 
filling of ILL, the v = 3 state can also be interpreted as 
a combination of = 2 and v = A states. 

Increasing the magnetic field beyond v = 5/2 forces the 
spin-polarized electrons in ILL to fall one- by-one down 
to OLL, starting from the electrons with highest angu- 
lar momenta. ILL maintains the spin polarization but 
the electrons in OLL tend to prefer spin compensated 
configuration and, as a consequence, there is a partial 
relapse towards lower total spin polarization. However, 
near v = 2 OLL gets a small spin polarization in high 
electron numbers, e.g., S' = 1 for = 48 and S' = 2 for 
A^ = 60 [see Fig. ^a)]. This result is in accord with the 
results of Ciorga and coworkers who found the collapse 
of the spin singlet state in large QD's at = 2 (Ref. Isl). 

The results at i/ = 5/2 can be generalized to larger 
particle numbers by assuming a spin compensated OLL 
and polarized ILL. The formula of i'oll (see Ref. 0) gives 
S « A^/10 for the spin polarization at i/ = 5/2, which is 
in a reasonable agreement with the calculated maximum 
spin polarization of the two largest dots in Fig. ^a). As 
seen above in the case of z^ = 7/2, we find also evidence 
of analogous polarization of the highest occupied LL in 
lower magnetic fields. This suggests that the results can 
be generalized to higher half-integer filling factor states 
(i/ = 9/2 etc.) in sufficiently large QD's. 

The predicted spin polarization of the states between 
2 < < 3 is interesting with a viewpoint on the spin- 
tronics in which electron spin degree is exploited to create 
novel device concepts and new functionality for electronic 
circuits. It is assumed that the alternation of states in the 
2 < ly < 3 regime leads to the familiar checkerboard pat- 
tern of conductance peak he ight s in the electron tunnel- 
ing through a QD device |lOllll| . However, the details of 
the current patterns are not completely understood. The 
checkerboard behavior is usually explained assuming os- 
cillations between the spin singlet S' = and triplet S = 1 
states jl^ which can be qualitatively modeled within the 
constant-interaction picture |l3l | with a phenomenologi- 
cal exchange term. Our SDFT calculations indicate that 
this model is sufficiently accurate in low electron numbers 
[see, e.g., data for = 12 in Fig. Gla)], but there are 
significant deviations from it for TV > 20. The constant- 
interaction model cannot take into account the correla- 
tion effects which give rise to the high spin polarization 
near the ly = 5/2 state in large QD's. This might affect 
electron tunneling probability through the dot for high 
electron numbers and, as a consequence, leave a trace to 
the pattern of conductance peak heights. However, the 
high probability transmission in the experiments occurs 
mainly through the edge region which is not spin polar- 
ized yi]j. The transport through the center of the dot 
gives rise to lower conductance peak heights, and evolu- 
tion of the ground state spin might leave a systematic 
pattern in the low probability transmission peak heights. 

We turn now into analysis of the exact many-body 
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wave function of the = 1/2 state. We use the ED 
method on OLL, and the results are compared to two 
models, the PF and the CF wave functions The PF 
wave function 'll is defined as 




with angular momentum L = N{N— l) — N/2. The value 
of the filling factor at < 1, i.e., beyond the maximum- 
density-droplet (MDD) domain is widely approximated 
to be = Lmdb/L. This definition is based on Jastrow 
states but the PF wave function has lyj which is slightly 
higher than 1/2. As the number of particles in the system 
increases, i^j for the PF approaches 1/2 from above as 
i^j = l/2 + l/2N+0{l/N)'^. The PF state is constructed 
for even TV, and we are computationally limited to TV < 8. 
The CF wave functions lal used here are defined as 




where Vfo is a determinant formed from the single-particle 
states ipn,i oc 2;""*"' 9", with I the angular momentum and 
n = 0, 1, 2 . . . is the LL index. The derivative is with 
respect to z, and the resulting many-body state has a 
total angular momentum L = N{N — 1) -|- X^il^i ^i- 

There are several possibilities for the single-particle oc- 
cupations in for the states near v,^ w 1/2. We con- 
sider two different CF wave functions: CFl that have 
the same angular momentum as PF states, and CF2 that 
have i/j = 1/2, so that L = N{N-1) = 2Lmdd- The first 
type have single-particle quantum numbers (n, /) of the 
occupied states as (0, i — 1) and (i, —i), i = 1, . . . , N/2, 
and in the second one the asymmetry in I quantum num- 
bers is recovered by moving a particle from {N/2, —N/2) 
to (1,0). The CF2 wave function is constructed also for 
odd A''. The numerical projection to OLL (derivatives in 
*o) is doable up to = 6. From these wave functions, 
CFl (and thus also PF) has an angular momentum that 
corresponds to the ground-state of A^ = 4, and CF2 for 
N = 5 and 6. In Table U we show the accuracy of the 
three wave functions discussed above. PF is clearly less 
accurate than CF which has very high overlaps. The rea- 
sonably high overlap (around .6) means, however, that 
PF cannot be an accurate excited state. Therefore we 
find no evidence of electron pairing in the 1^—1/2 state. 

Figure 121 shows the conditional wave functions j^, ob- 
tained by fixing A^— 1 electrons and moving the remaining 
electron, for the PF and corresponding ED state. One 
can see that the most probable positions are different 
in ED and PF and the ED conditional density is more 
localized than the PF one. 

In the 2DEG the electrons in the half-filled ILL can be 
described by a = 1/2 state on OLL having a screened 



TABLE L Overlaps of the trial wave functions with exact one 
for different particle number TV. 



_/V PF CFl CF2 

4 0.922124 0.999903 0.997978 

5 - - 0.998921 

6 0.789996 0.993461 0.996096 
8 0.586370 




FIG. 2: (color online) Electron positions (crosses), conditional 
electron densities (contours), and phases (gray-scale) of the 
N = 8, L = 52 state from ED (a) and PF model (b and c). 
We probe with the rightmost electron and densities are on 
logarithmic scale. The phase changes from n to — vr on the 
lines where the shadowing changes from the darkest gray to 
the white. In (a) and (c) electrons are on the most probabe 
positions, and in (b) we use ED coordinates for PF. 

electron-electron interaction. This state is accurately 
modeled using a PF wave function P, |^. The mean- 
field results above suggest that because of the different 
structure of the 1/ — 5/2 states in QD's and the 2DEG 
systems, the v = h/2 state is not related to v — 1/2 in 
parabolically confined QD's. The ILL in QD's is com- 
pact and does not have the required angular momentum 
of the V = \/2 state. One might assume that realistic 
inter-electron potentials, which are modified by screen- 
ing from the electrodes and finite thickness of the sample, 
would induce electron pairing as observed in numerical 
studies of the v — h/2 state in the 2DEG. However, ED 
calculations with a screened electron-electron interaction 
indicate that overlaps of the PF with the exact state in- 
creases only marginally. Therefore we find no evidence of 
electron pairing in the analyzed half-integer filling factor 
states in QD's confined by a parabolic external potential. 

The confinement potential in actual QD's may vary 
considerably from the parabolic shape [1J|. However, 
this fact seem not to affect notably the appearance of the 
V — 5/2 state in QD's. Our SDFT calculations for a 60- 
electron QD defined by an infinite potential well of radius 
150 nm yield the v — h/2 state similarly to the parabolic 
case. As the main difference, the state is not separated 
into two V ~2 and = 3 domains in the infinite- well QD. 
Instead, electrons on OLL are strongly localized near the 
edge, whereas the spin-up electrons on ILL contribute 
in the mid-region leaving the electron density in the core 
relatively flat for both spin types. In the ground-state en- 
ergetics the behavior between different QD's is similar. 
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FIG. 3: (color online) Ground-state magnetization (upper 
curves) and total spins (lower curves) of 60-electron quantum 
dots in parabolic (a) and infinite-well (b) potentials. The fi- 
nite size counterparts of the integer and fractional quantum 
Hall states are marked in the figures (dashed lines). 

as we demonstrate in Fig. 13 showing the SDFT result 
of the magnetization, M — —dEtot/dB, for 60-electron 
parabolic and infinite-well QD's, respectively. The peaks 
in M coincide with the integer filling factors and consti- 
tute the finite size counterpart of the de Haas-van Alphen 
effect. In both cases, the v = 5/2 state leaves cusps in 
the magnetization, and the data shows also possible de- 
velopment of other half-integer filling-factor states such 
as = 7/2 in lower magnetic fields. The results are con- 
sistent with the few direct magnetization measurements 
of ensembles of large QD's which have shown a cusp-like 
structure near 1^ — 5/2 (Kei. Il^ll . Visible signatures are 
also expected to be found in the chemical potentials of 
electron tunneling experiments. 

To conclude, we predict emergence of states in quan- 
tum dots which can be regarded as finite size counter- 
parts of the half-integer filling factor states in the 2D 
electron gas. The highest occupied Landau level of these 
states is found to be spin-polarized and the formation 
of them leaves signatures in the ground state energetics. 
The numerical results indicate that these states share 
many of the characteristics of their infinite size counter- 
parts in the 2D electron gas, although we find also major 



differences. Due to observed stability of these states we 
postulate that the formation of them is a general prop- 
erty of finite fermion systems in external confinement. 
In parabolic confining potential the composite fermion 
picture is found to be a more accurate description of the 
= 1/2 state than the Pfaffian wave function, which sug- 
gests that electrons are not paired in this geometry. The 
electron pairing, which is predicted to occur in the two- 
dimensional electron gas, might be recovered in quantum 
dots in the limit of very large electron numbers and with 
a weak non-parabolic external confinement potential. 
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